Wave scattering by small particles in a medium. 
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Abstract 

DDD 

Wave scattering is considered in a medium in which many small particles are 
embedded. Equations for the effective field in the medium are derived when the 
' number of particles tends to infinity. 

o ■ 

^ ■ 1 Introduction 
o 

^ ! Assume that waves in the medium are described by the equation 

r-| ■ Lu := [{aip{x)u^p)i + k'^n{x)]u = in M^, 



> 



where over the repeated indices summation is understood, u^p := and the Green 
function, satisfying the radiation condition, solves the equation: 



^ I LG = —S{x — y) in 



If there are M particles D^, placed in the medium, situated in a bounded domain D, 
outside of which 

aip{x) = 5ip, n{x) = 1, x e D' ■.= W^\ D, (3) 

where ajp(x), n{x), are C^-smooth functions, dip is the Kronecker symbol, and the ellip- 
ticity condition holds: 

3 3 3 

Cl I^PI^ - 5Z ^iP^P^i ^ ^2 > 0, 

p=l c,p=l p=l 
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where t G is an arbitrary vector, then the scattering problem consists of solving the 
equation 



M 



LU = in R3\|J Dm, (4) 



m=l 



U\g = on Sm, l<m<M, (5) 

M 

U = Uq + Y^ / G{x,s)am{s)ds. (6) 

m=l 

Here LUq = 0, uq is the scattering solution in the absence of particles, i.e., if M = 0. 

By Ramm's lemma |3j, p. 257, one can define the scattering solution Uq in the absence 
of small particles by the relation: 

G{x,y) = g{y)Uo{x,m + o{l)], \y\ ^ oc, ^ = "A (7) 

where g{y) := We assume that ka << 1, where a = |maxi<m<M diamD^- 

The aim of this paper is to develop a general approach to wave scattering in a medium 
in which many small particles are embedded. Smallness of the particles is understood in 
the sense ka << 1. The functions n{x) and assumed practically constant on 

the scale of the wavelength 

k{\Vn\ + iVflipl) << 1. 

We generalize the approach developed in [5], [H], [Z]. Earlier works are [T], [TU], [H], 
to mention a few. 

Our basic result is a formula for the wave field in the medium in which small particles 
are embedded. This field solves equation (jl])-© and satisfies the radiation condition at 
infinity: 

^^^'^'^ -tk{U-Uo) = o(-), r:=|xHoo. (8) 



We assume that 

(i>a, /ca<l, (9) 

where d = min^^j dist {Dm,Dj), and the dist denotes the distance between two sets. 
Near any point x E D, such that mini<m<A/ dist {x, Dm) » 1, one calculates the wave- 
field U by the formula 

M 

U{x) =Uo{x) + '^G{x,Xm) [I + ikiy ■ {s - Xm)]crrnis)ds, (10) 

m=l ■^'^™' 
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where ^ \x-^\ ' ^^"^ ^ depends and on the functions n{x), aip{x), The term 

iku ■ {s — Xm) comes from the formulas 

G{x, s)am{s)ds = G{x, Xm) { I (Tm{s)ds + I - — ^ ' J^^^ ^ ' — '^^ds 

iku ■ (s — Xr 



G{X, S) - G{X, Xm) _ Jo ^yG{x, Xm + r{s - Xm)) " (s - Xm)dT 
G(^Xy Xm) ^(*^; •^m) 



where ly = v{x, Xm), ^ = jf^^^Efj • may consider z/ as a known vector because G{x, y) 
is known. 

In a generic case, when Qm '■= fo c^m{s)ds ^ 0, the assumption ka « 1 allows one 
to neglect the term ikv ■ s = 0{ka) in (flOj) and to write (ITOl) as 

M „ 

U{x) =Uq{x) + ^^G{x,Xm)Qm: Qm = CTmds, (11) 

m=l 

If Qm = 0, then the term iku ■ s cannot be neglected. We discuss this case in Section [31 
A physical example of such a case is the scattering by acoustically hard particles when 
the boundary condition is the Neumann one: =0, l<m< M. 

2 General methodology 

Let us first assume that ( fTTj) is applicable and calculate Qm- In a neighborhood of Sj 
one has the exact boundary condition ([5]), which can be written as: 



f G{s,t)aj{t)dt = - {uo{s) + ^G{s,Xm)Qm\ 



Ue{s), (12) 



where s G Sj and Ue is the effective field acting on Dj. The basic assumption is: 
We assume that Ue{s) is practically constant on the distances of order a. 
As |s — tl -^0, one has 

G{s, t) = ^^^^ [1 + 0{ka)], \s-t\^0, (13) 

where s,t G Sj and we have assumed for simplicity that Oij = 6ij. In the general case one 
replaces the function go{s,t) := ^^|^_^| on the surface Sj by the fundamental solution of 

the operator Ylip=i ^ipi^j) dx dx ' which can be written explicitly and analytically: 
G{x,y) 



4n./d^tM [a'i^'\x,-y,){xp-Up)Y/^ 
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where the matrix a^^"^^ is inverse of the matrix aip{x), Xi is the z-th Cartesian component 
of the vector x (not to be confused with the vector Xj = x). 
If Qip = 6ip, then G solves the equation 

G{x,y) = go{x,y) + I g{x,t)n{t)G{t,y)dt, go{x,y) = -. (14) 

Jd 4:7r\x-y\ 

The integral in (|T4|) is 0(1) as ja: — y| ^0, so f|T3l) follows from f|T^ . We have: 

G{x,y) = go{x,y)[l + 0{\x-y\)], \x - y\ ^ 0. (15) 

In this paper > is assumed fixed. The error term in f|T3|) is 0{ka), and the small 
parameter ka is dimensionless. Replacing G{s, t) in f|T2l) by its expression f|T3|) and 
neglecting the small term 0{ka), one gets the following integral equation for aj: 



Si 



go{s,t)aj{t)dt = -Ue{xj), (16) 



where we may replace Ue{s), s E Sj, by the quantity Ue{xj) because |xj — s| < a, and 
Ue was assumed practically constant on the distance of order 0(a). Equation flTBl) is an 
equation for the electrostatic charge density crj{t) on the surface of a perfect conductor 
Dj, charged to a constant potential —Ue{xj). The total charge on the surface of this 
conductor is 

Qj = / ajdt. 
Js, 

One knows from electrostatics that the total charge on the surface of a perfect conductor 
Dj equals to the product of the electrical capacitance Gj of this conductor and the 
potential —Uf.{xj), to which this conductor is charged: 

Qj = -CjUe{Xj). (17) 

In ^ p. 26, formula 5.12], analytical formulas are derived for calculating the electrical 
capacitances of conductors of arbitrary shapes with any desired accuracy. 

Thus, one may consider the capacitance Gj to be known. Equation flTTl) can be written 

as 

M 

U{x) =Uq{x) - y^^G{x,Xm)GmU{xm), dist (x, D„) > o. (18) 

Equation (fT7|) can be considered as a linear algebraic system for finding Q^. Namely, 
^ and (dZ]) imply: 



Qj = ~Gj I ^o{^j) + ^^^i' ^rn)Qm j 



(19) 
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Equation fll9l) is a linear algebraic system for the M unknowns Qj. The matrix elements 
of this sytem are known because the Green's function G{x, y) is known. If the condition 



C,Y,\G{xv^m)\<l, (20) 

holds, then system f|T9|) can be solved by iterations and the iterative process 

^ _c^Uo{x,) - Y,Cfi{x„xMt\ Qt^ ■■= -CMx,), (21) 

converges to Qj at the rate of a geometric series as n oo. If M — ^ cxd and the limit 

/ C{y)dy = lim J] C, (22) 

Jd M^oo 

D,CD 

exists for any subdomain D G D, where the function C{y) is integrable, then this function 
is the limiting density of the electrical capacitance of the small particles. If the limit (122!) 
exists, then equation (fTSll in this limit takes the form: 

U{x)=Uo{x)- [ G{x,y)C[y)U{y)dy. (23) 
Jd 

Here U{x) is the limit of Ue{x) as M — > cxd and it is assumed that (1221) holds. 

Applying the operator V + k'^n{x) to ( l23ll and using equation ([21) yields the following 
equation: 

[V^ + Pn(x)]W = Cix)U{x). (24) 
This is a Schrodinger-type equation for the effective (self-consistent) wave field U: 

[V^ + k^~ q{x)]U = 0, (25) 

where 

q{x) ■.= C{x) + k'^[l-n{x)]. (26) 

Thus, in the limit M ^ oo under the assumption fl22]) . we have derived a linear 
Schrodinger-type equation (125!) for the effective wave field U. If the small particles are 
identical, the electrical capacitance of a single conductor with the shape of a particle is 
Co, and the density of the number of small particles is N{x), then C{x) = N{x)Co, so 
N{x) can be calculated if C{x) and Cq are known. By the density N{x) of the number 
of small particles we mean the function N{x), defined by the formula, similar to 



/ N(x)dx = lim > 1. 
Jd m^'^ 



D.CD 
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Let us summarize our method for calculating the effective field in the medium in 
which M small particles are embedded and the Dirichlet condition holds on their 
surfaces. 

This field is calculated by formula ffTTj) at any point x such that min +m dist {x, Dm) >> 
a, and the unknown numbers Qm are calculated by solving linear algebraic system f|T9|) . 
If condition fl2(l holds, then system (HM is uniquely solvable by iterations. Condition 
f l20p always holds if M is fixed and a is sufficiently small, because Cm = 0{a). 

If M — s> oo and condition (122!) holds, then the field U can be found from the integral 
equation (!23|) . Solving this equation is equivalent to solving equation (!25l) where U —Uq 
satisfies ([8]). The function n{x) is known, therefore G{x,y) is also known. 



3 The case when Qm = 0. 

Consider the same problem as in Section [2] but with the Neumann boundary condition 



dU 

m 



0, 1 < m < M, 



(27) 



in place of the Dirichlet one. Here is the unit exterior (i. e. , pointing outside of Dm) 
normal to Sm, 1 < m < M. We look for the solution of the same form ([6]), and the 
methodology is the same, but now Qm = 0, 1 < m < M, as we show below. Thus, the 
term iku ■ s in (fTOj) becomes important. Arguing as in Section [21 we obtain in place of 
( IT2|) the following equation 



Ou ( s) 

a,{s)=Ajaj + 2^^, s e S„ 



where 



A,a, := 2 [ 



dN ' 
dG{s,t) 



and the known formula from the potential theory: 

d 



ON: 



G{x,t)aj{t)dt 



(Jj{t) dt, 



(28) 
(29) 

(30) 



was used. Here is the normal derivative on Sj from the outside of Dj. 



If \s-t\ ^0, then 



G{s,t) ^ go{s,t) 



1 



47r\s - t\ 

where we again assume for simplicity that aip = 6ip. We have 

G{Xj^ s) G{Xj^ '^m) 



J , ^) (—T^Xj , X 



1 + 



s e S„ 



(31) 
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Since \s — Xm\ < a and \xm — Xj\ ^ a, we have 



Thus, 

(jT^^XjjSJ Kj\Xj^Xm) -If \ fnn\ 

^ ^ikv ■ {S- Xra), (32) 

(j \Xj , X^ j 

where, by the mean value theorem, one has 

iku = Vyg(a:,^m + r(g-x^)) ^ < ? < 1. (33) 

If \xj — Xm\ » a, then 

I ^ I I J I 

and depends on the functions n and Cjp. Since n(x) is known, then G{x, y) is known, 
so ^{xm) can be considered as known. 

If ttip = 6ip, then the integral equation for G can be written as 

G{x, y) = gix, y) - ! g{x, Oq{i)G{i, y)d^ := ^ - TG, (35) 

J D 

where T is the integral operator, defined in (!35|) . 

^ik\x—y\ 

g{x,y):=— q{y) := k^[l - n{y)], q = in D' ■.= R^\D. 

A7r\x — y\ 

The operator T is compact in L'^{D). 
The function 

_ -f 7 1 1^1 

^ z/c-j if fc ?/ — 1. 

fi' \y-x\ 

More precisely, we have used the exact formula: 



^y9 = 9{x, y) {ik - ^ ) ■ ^ ^ 



\y — x\ \y — x\ 

and have neglected the term -p^^^ compared with \ik\ = k. 
This is justified if k\x — y\ ^ 1. 

Equation fl35|) is uniquely solvable in L'^{D) because it is of Fredholm type and its 
homogeneous version has only the trivial solution. Indeed, if h solves the homogeneous 
equation (1551) . then 

{V^ + k'^ - q{x))h = inM^ 
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where the potential q is compactly supported, A;^ > 0, and h satisfies the radiation 
condition. It is known that this implies h = 0. In Kato's paper [2J a similar but much 
stronger result is obtained: it is not assumed that q is compactly supported, the potential 
may staisfy the assumption \q{x) \ = o{\x\~^) as |x| —>■ oo. 
The unique solution to (|35|1 has the form 

G={I + T^)g, 

where / is the identity operator and Ti := (/ + T)^^ — / is a linear compact operator in 
L\D). 

Therefore 

VyG ^ (j + rpv,^ 

G G{x,y) ■ 

li\x — y\ ^ oo and A; > is fixed, then 

^^""'^^ -l + o(—^—^ dist(x,D)>l, dist(i/,D) > 1. (36) 



g{x,y) \\x-y 
This follows from equation (l35ll if one takes into account the following estimates: 



\g{x,^)\ = 0[ i^^^i ) as |x-^|^oo, 
1 



\G{i.y)\ = 0{j-^A as le-l/l-oo, 



and 



J-=l-. 1^ \ = 0{y^X \x-y\^oo. (37) 

Jd f - A\z -y\ \\x-y\J 

The last estimate holds if D C is a bounded domain. To prove estimate fl57|) . take 
the origin at the point note that (1371) holds if dist (?/, D) < diamD, and J decays 
when dist (y, D) grows remaining less than , if the gravity center of D moves along 
the line joining x and y. 

If X or y remain in D and \x — y\ oo, then the relation ^^^'^1 ^ 1 does not hold, in 

f I »i 5 g(x,y) ' 

general. 

Since G{x,y) is known one may assume that = \iy{x)\ = \iy{x,y)\ is known. 
Because of ( IT3|) , the operator Aj in equation (!28|) can be approximated by the operator 

It is known (see ^ p. 96, formula (7.21)]) that 

Aads = - / ads. (38) 



3 "J 
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Therefore equation (125]) implies 



f ajds = I -^ds = f Auedx ^ VjAue{xj), Vj := \Dj\, (39) 

J Sj J Sj J Dj 

where \Dj\ is the volume of Dj. We had assumed that Ue{x) is C^-smooth. The term 
Vj = 0{a^), and Aue = 0{k^), because AWq = -Pn{x)Uo = 0{P). Thus, 

VjAue = 0(A;V). 

Let us show that the second term in (ITU]) , namely, ik ly ■ {s — Xj)ads, is of the same 

order of magnitude 0(/c^a^). 

We have assumed that z/ is practically constant on the scale of order 0(a). Thus, 

^ r ^ ^ 

ik'^Ug [s - Xj)qaj{s)ds = -ik ^o^jl^qp-^- (40) 



l,p=l 

Here we have used the following result (see [H, p.98]): 
If 

hp = Ahp - 2Np, 

then 

Xj)qhpds Vjf3pq , 

where (s — Xj)q is the q—ih Cartesian component of the vector s — xj, Vj is the volume 
of the domain Dj, and Ppg is the magnetic polarizability tensor, defined in [5] p. 62, 
formula (5.62)]. In the cited formula one takes the magnetic constant /io = 1, and the 
origin is at the point Xj, the gravity center of Dj. The right-hand side of (1^0]) is of order 
O(a^) = 0{k'^a^) because > is fixed. Thus, if the Neumann boundary condition holds 
on Sm, 1 < m < M, then the terms a^ds and ik v ■ {s — Xm)ds in (ITU]) are of the 
same order of smallness as ka — >■ 0, and > is assumed fixed. 

Let us compare this conclusion with the case when the Dirichlet boundary condition 
holds on Sm, 1 < m < M. In this follows from (|T7]) . Qm = 0{a), because 

Cm = 0{a). If a ^ 0, then 0(a) > 0{a^). This is the reason for the different physical 
conclusions in two cases. In the case of the Dirichlet boundary condition the scattering 
amplitude is of order 0(a) and the scattering by a single small particle is isotropic, while 
in the case of the Neumann boundary condition the scattering amplitude is of order 
0(A;^a^) and the scattering by a single small particle is anisotropic. 

Let M — i> oo in the case of the Neumann boundary condition, and assume a^p = Sip. 
Equations ([ID]), ([211) and (HOD yield 



U{x)=Uo{x)+ / G{x,y) 



D 



^ dU 



AU{y) -ikY^ Ug{x,y)—(3pg{y) 



p,q=l 



ip 



v{y)dy. (41) 
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Here we have assumed that for any sub domain D G D one has: 

/ v{y)dy = hm ^i' (42) 



D.CD 



/ Ppgiy)dy= lim V 
Jd m^oo 

DjCD 

and then we have passed to the hmit M — cx) in equation (ITOl) . 

In in p. 55, formula (5.15)], analytical formulas are given for calculating tensor /3j 
with any desired accuracy for a body Dj of arbitrary shape. 



(i) 

■pq 



4 Conclusion. 

The methodology we have developed for solving many-body wave scattering problem for 
small particles, embedded in a known medium, has the following new features: 

1) Scattering by small particles is considered in a medium. 

2) If the number M of the particles is not very large, the scattering problem is 
reduced to solving linear algebraic systems with matrices, whose entries have physical 
meaning. Analytical formulas for calculating these entries are obtained in [3]. The 
reduction to linear algebraic systems bypasses any usage of integral equations, which are 
usually serving as a basic tool in the scattering theory. 

3) If M oo, then some integral equations (equations fl23|) and fHTj) ) are derived for 
the effective field in the medium in which small particles are embedded. 

4) In the case of the Dirichlet boundary condition the relative volume of the embedded 
particles tends to zero as a — 0, M — cx). 

Indeed, the number of small particles per unit volume of the medium is O(^), where 
d is defined in (I9]), the volume of a single particle is O(a^), so the relative volume of the 
small particles is O(^) ^ as M — oo, since a — >• if M — oo. 

In the case of the Neumann boundary condition the relative volume of the small 
particles tends to a finite non-zero limit v{y) (cf. fl42p ). 

5) The methodology, developed in this paper, can be used in the problems of electro- 
magnetic wave scattering by small particles embedded in a known medium. It can also 
be used in some nanotechnological problems, consisting of creating "smart" materials 
with the desired properties, for example, wave- focusing properties (see [8j, ^). 
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